



III ' lllfll W$i |S!::i 

g§ V'.; M :.< Artlmr:.!,.- GlassmaB.' li ii Mig ;■' ||gj| g ; = . aA;-' ■■H^^ r: |||p|| |||| Jg ||l|i||g| 

1 . ■-' GRANT NA03-1165- \ ' ' I ' - ' - 

. ; | | . FEBRIJ ARY I.9&2 | , | . . ' |j ' •/ | | V | | Jj g| | 




■ r ' : V"': '■:;■;?■■ v'.l' V.A ; :i v ^: \:' /Vu" i : ■ - :■ " r ;" ; V : .'.'•& 1 -s -i 1 -'; ■ ■■':■ -"■". /''t: . ■ :i< 


(MASA-CR-4430) COMPUTER CODE FOR 
PRELIMINARY SIZING ANALYSIS OF axial-flow 
TURBINES Final Report CToleao Uni v» ) 


N9 Z - 20 1 9 o 



NASA Contractor Report 4430 


Computer Code for Preliminary Sizing 
Analysis of Axial-Flow Turbines 


Arthur J. Glassman 
University of Toledo 
Toledo, Ohio 


Prepared for 

Lewis Research Center 

under Grant NAG3-1165 


l\J/\SA 

National Aeronautics and 
Space Administration 

Office of Management 

Scientific and Technical 
information Program 

1992 



Summary 

This report presents a computer program for the preliminary 
sizing analysis of axial-flow turbines. The computations are 
based on mean-diameter flow properties and a stage-average 
velocity diagram. Given as input to the program are power 
or pressure ratio, mass flow rate, inlet temperature and 
pressure, rotative speed, inlet and exit diameters (either hub, 
mean, or tip), exit radius ratio or stator exit angle, turbine 
loss coefficient, and gas properties. Computations are then 
performed for any specified number of stages and for any of 
three types of velocity diagrams (symmetrical, zero exit swirl, 
or impulse) or for any specified stage swirl split. Exit turning 
vanes can be included in the design. The program output 
includes inlet and exit annulus dimensions, exit temperature 
and pressure, total and static efficiencies, flow angles, and 
last-stage absolute and relative Mach numbers. 

The analysis and code presented herein are modifications 
of those in reference 1. New features added to improve 
modeling rigor and extend code applicability include a gener- 
alized velocity diagram, a more flexible meanline path, a reheat 
model, a radial component of velocity, and a computation of 
ffee-vortex hub and tip velocity diagrams. Also, a loss-coefficient 
calibration was performed to provide recommended values for 
current technology airbreathing engine turbines. 

The analysis method and a complete description of input and 
output are presented in this report. Sample cases are included 
to illustrate the use of the program. 


Introduction 

The preliminary analysis of a power or propulsion system 
involves many repetitive calculations to determine system 
performance, component performance, and component geom- 
etries over a range of conditions. This must be done in order 
to eventually determine the best system and operating conditions. 
For this type of screening analysis, complete design accuracy 
and detail for the components are not necessary. Approximate 
and rapid generalized procedures, rather than complex and 
time-consuming detailed design procedures, are sufficient to 
yield the desired component overall geometry and performance 
characteristics. 

A computer program for the preliminary sizing analysis of 
axial-flow turbines was presented in reference 1. This meanline 
analysis uses the efficiency computation procedure of refer- 


ence 2, which is based on a stage-average velocity diagram, and 
incorporates assumptions consistent with turbine technology 
from about 1970, when the program was developed. With 
turbines for some current and future propulsion systems 
exhibiting higher pressure ratios and more highly inclined 
flowpaths than in the past, it became apparent that modifica- 
tions were required to improve modeling rigor and extend code 
applicability. New features added to accomplish this include 
a generalized velocity diagram, a more flexible meanline 
flowpath, a radial component of velocity, and the computation 
of free-vortex hub and tip velocity diagrams. Also, a loss- 
coefficient calibration was performed to provide recommended 
values for current technology airbreathing engine turbines. 

The analysis method and a complete description of input and 
output are presented in this report. Sample cases are included 
to illustrate use of the program. Because reference 1 is no 
longer readily available, this report describes the entire 
analysis, not only the modifications. 

Symbols 

A stage loss parameter 

A an annulus area, m 2 ; ft 2 

C blade loss parameter 

C A dimensional constant, In rad/rev; 60 sec/min 

C B dimensional constant, 1; 550 ft • lb/sec • hp 

c x axial chord, m; ft 

c p heat capacity, joules/kg ■ K; Btu/lb * °R 
D diameter, m; ft 

E squared ratio of stage-exit axial velocity to stage- 
average axial velocity 

F blade loss weighting factor 

g dimensional constant, 1; 32.2 ft • lbm/sec 2 * lbf 
Ah specific work, J/kg; Btu/lb 

i stage number, i = 1,2,. . .,n 
J dimensional constant, 1; 778 ft • lb/Btu 
j stage number for change in meanline slope 

K turbine loss coefficient 

k coefficients for axial chord correlation 

M Mach number 
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N rotative speed, rad/sec; rpm 

n number of stages 

P shaft power, W; hp 

p pressure, N/m 2 ; lb/ft 2 

R gas constant, J/kg K; ft • ibf/lbm °R 
Re Reynolds number 

r radius, m; ft 

T temperature, K; °R 

U blade speed, m/sec; ft/sec 

V absolute gas velocity, m/sec; ft/sec 

W relative gas velocity, m/sec; ft/sec 
w mass flow rate, kg/sec; lb/sec 
a absolute-flow angle from axial direction, deg 
# relative-flow angle from axial direction, deg 
7 specific heat ratio 

rj static efficiency 

7 }' total efficiency 

6 meanline angle from horizontal, deg 
p viscosity, N • sec/m 2 ; lb/sec • ft 
p density, kg/m 3 ; lb/ft 3 

\J/ work factor 


Subscripts: 

a first stage 

abs absolute 

cr critical 

ev exit vane 

ex turbine exit 

h hub 

i stage i, i = 1,2,. . . 

id ideal 

in turbine inlet 

is value at turbine inlet entropy 

m mean section 

me meridional component 

n last stage 

rel relative 

ro rotor 

st stator 

t tip 

u tangential component 

x axial component 


0 stator inlet 

1 stator exit 

2 rotor exit 

Superscripts: 

~ turbine overall 
7 absolute total condition 

" relative total condition 


Method of Analysis 

The method is based upon an analysis of the flow at the 
turbine mean diameter. Radial gradients of the flow properties 
are not included in the flow and efficiency analyses. However, 
free-vortex hub and tip velocity-diagram parameters are 
computed in order to. indicate the possible severity of flow 
conditions at the endwalls. The specific heat ratio is assumed 
constant throughout the turbine. For any given turbine, all 
stages, except the first, are specified to have the same shape 
velocity diagram. The first stage differs only in that the inlet 
flow can he at any specified flow angle. The velocity diagram 
shape depends upon the stage work factor value and the 
specified velocity diagram. Three specific types of velocity 
diagram (symmetrical, zero exit swirl, and impulse) or a 
diagram with a specified swirl split (i.e. , distribution of swirl 
velocity between rotor inlet and rotor exit) can be used. 
Figure 1 illustrates the velocity diagram symbolism and 
presents the three specific types of velocity diagram for three 
values of stage work factor. The fundamentals of velocity 
diagrams and of the flow and loss modeling used herein can 
be found in reference 3 . 

Various input options dictate the exact nature of the 
calculation procedure. There is, however, one basic procedure 
that is direct and without iteration. This basic procedure will 
be presented, and then the alternative procedures required for 
the various input options will be discussed. 

Basic Calculation Procedure 

The required inputs for the basic procedure are shaft power 
P, mass flow rate vv, inlet total temperature T' n , inlet total 
pressure p [ n9 rotative speed V, inlet mean diameter D m in> exit 
mean diameter D m ex , stator exit angle a h gas constant R, 
specific heat ratio 7, viscosity /a, loss coefficient K, and 
the squared ratio of stage-exit to stage-average (i.e., {V me Q 
H- V me l + V me ,2 )/ 3 ) meridional velocities E. Also specified 
for each calculation are the number of stages n and the type 
of velocity diagram. For a multistage turbine, the input variable 
specified as inlet diameter is used to calculate first rotor blade 
speed and annulus dimensions at the first rotor exit. Therefore, 
it is truly an inlet diameter only if the hub and tip diameters 
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(a) 



(a) Symbolism. 

(b) Effect of diagram type and work factor on velocity diagram shape. 
Figure 1.— Velocity diagrams (constant axial velocity). 


are assumed constant across the first stage. For a one-stage 
turbine, the specified exit diameter is used for the calculations, 
and the inlet diameter is of no significance. 

The first- and last-stage blade speeds are 


T r * ND mAn 

a r 

U A 

(i) 

jj _ ^ ND m,ex 
" Ca 

(2) 


For more than two stages, there are three options for stage 
blade-speed variation (i.e., mean flowpath): 

(1) Linear variation between first and last stages 


Ui = — — ^ (i -1) + U a (3a) 

n — 1 


(2) Constant from first stage to j th stage (i = 1 to j) 
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Ui = U a 


(3b- 1) 


and then linear to last stage ( i = j to n) 


U, = — — (i-j) + U a (3b-2) 

« -/ 


(3) Linear from first stage to j th stage (i = 1 to j) 

Ui = — f a a - 1 ) + U a (3c- 1) 

J - 1 

and then constant to last stage (i = j to n) 

Ui = U n (3.C-2) 

Turbine specific work is 


— CnP 

A h'=~f~ (4) 

J w 


and is equal to the sum of the specific work of the stages: 

( 5 ) 

Expressing stage specific work in terms of stage work factor as 

( 6 ) 

and substituting equation (6) into equation (5) yield 
V Uj 

Ah' =L —i (7) 

i=l gJ 


Ah' = E Ah- 

i=l 


Uj 

A h!=^J 

gJ 


£ 

I 

& 



pV+i 


Entropy 

Figure 2. -Temperature-entropy diagram for * th stage. 

additional features: (1) the turbines considered in this report 
are not restricted to a constant mean-section diameter, (2) exit 
vanes to provide axial flow leaving the turbine can be included 
in the design, (3) the velocity diagrams also can be specified 
as symmetrical or by an input swirl split, and (4) the effect 
of turbine reheat is included. 

Referring to figure 2, turbine overall total efficiency is 
defined as 


Ah' 


Ah' 


Ah' 


n ~\u^ n " /T' \ 

** Em,-, 

1 1 V T U 


(9a) 


or alternately expressed as 


Since the velocity diagram shape is specified to be the same ^ 

for all stages, the stage work factor is the same for all stages q ' = ■ — — 

and is computed by rearranging equation (7) as ^ / Ah- d \ / Ahh 

i KlfAwJKKi') 


xf, = 


gJ Ah' 

n 

£ u} 

/= 1 


( 8 ) 


The value of \p is the primary factor determining turbine 
efficiency. 

The method used for computing turbine efficiency is basically 
similar to that presented in reference 2, but has the following 


Combining equations (6) and (7) yields 
Ah; uf 



1 


and, by definition 


(9b) 


(10) 
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( 11 ) 


Vi 


Ah; 

Ahjd,i 


Substituting equations (10) and (11) into equation (9b) yields 


Some of the terms within the parentheses in equation (17) 
are the same for all cases, while others depend on stage 
location, velocity diagram type, and use of exit vanes. The 
rotor weighting factor F w and rotor loss parameter C w are the 
same for all cases: 


V =-? 



( 12 ) 


F rn = 2 


(19) 


C ro = 2 cot oq 1-^1 + I ^ - - ) + , 

\A Vj \AV U +J VAK„ f 


, in,iY ± /^, iv . (v u , 2 i 


The terms T{ and T is l , using equation ( 6 ), can be written as 

AhU 


77= T t . 


T t - , 


uU 


* 


Ah is,i— I 

'T't 'T'f __ T/ 

/ is,i * is,i — i 1 is, i — 1 


gJCp 

T U-‘\f 1 


03) 


* \ I I t 

T{-\ ;\m'-\)gJc p 


where 


R 


7-1 J 


(14) 


(15) 


Since both T( and T[ are equal to T{ n at the turbine inlet 
(i.e., at i = 1), equations (13) and (14) can be evaluated by 
recursion once the 77 / are known. 

The stage efficiency computation method is explained fully 
in reference 2 , and only the key equations are presented here. 
Stage total efficiency can be expressed as 


1 


A , 

1 + — t/' 
2 


(16) 


where the stage loss parameter A is 


A =■ 


KRq 


-0.2 


cot Ci\ 


(FstCst + F W Cro + Qv) 


(17) 


The constant of proportionality K , called the turbine loss 
coefficient in this report, must be determined empirically, and 
recommended values based on experimental efficiencies are 
presented later in this report. The Reynolds number used in 
this calculation is defined as 


Re = 


2w 




(18) 


( 20 ) 


For last stages of turbines having exit vanes, 


^ . , 2 ,'VuA 2 ,( V ».2 X2 

C ev = 2 cot a i ( | 4- 


A VJ \AV U ' 


( 21 ) 


For all other stages 


C ev = 0 


( 22 ) 


For first-stage stators having a specified inlet angle a 0 , the 
stator loss parameter is expressed 


F'St ~ 


1 + cot 2 «i (2 + tan 2 a 0 ) 




(23) 


For all other stators 


( 1+2 cotV,)( ^ ) + 


U, 2 

AK, 


(24) 


The stator weighting factor F st also depends on whether or 
not the stator is a first-stage stator. For a first-stage stator 
having a specified inlet angle a 0 , 


K,= 


3 tan a 0 
tan a, 

tan a 0 


(25) 


tan a 


For all other stators 


1 - 


F, = - 


3(F„, 2 /AFJ 

Vu.i/AV, 


(26) 
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The rotor inlet and exit swirl parameters V Utl /AV u and 
V u2 /AV u of equations (19) to (26) depend on the specific type 
of velocity diagram or are input for the general diagram. The 
following table presents the relations for evaluating the swirl 
parameters: 


'TV J>/ 

1 ex l in 


Ah' 


( 35 ) 


T = T' - 

* ex 1 ex 


IgJCp 


(36) 


■ 

Velocity 

Inlet swirl 
parameter, 

Ka'AK 

Exit swirl 
parameter, 
V« 2 IW U 


1 + 

1 - 1 

Symmetrical 



2i 

2<P 

Zero exit swirl 

1 

0 


2 + \J/ 

2 -lA 

Impulse 





General 

Input 

AV « 


With turbine total efficiency obtained from the foregoing 
equations, the turbine exit velocities and state conditions are 
computed as follows: 


AK„,„ = *U n (27) 

V u ,l,n ~ AV u ,n (28) 


Pex = Pex 



(37) 


In order that the turbine total and static efficiencies be 
consistent with the computed exit velocity, the static efficiency 
is computed as 


Ah ' 


V = 


c T 

L 'P* in 




(38) 


In order to calculate continuity for a flowpath having significant 
inclination, it is necessary to determine a flowpath slope. 
Consistent with the preliminary nature of this analysis, a 
correlation relating axial chord to mean diameter was used 
to determine axial length. The correlation, based on the geometry 
of existing turbines, is of the form 


V u .2 ,n = 


y.u.2 

av u 


AV„ 


(29) 


c x — k\D m 4- k 2 D 2 m + k 3 + k 4 D m 1 (39) 


V/ne.n l.nCOtO!] 

v = sTf v 

v me, ex vt/ r me,n 



no exit vanes 
exit vanes 


(30) and was established over a range of mean diameters from 0.10 
to 1.25 m (0.33 to 4.17 ft). Because of the scatter in the 
available data, correlations were determined for low, mean, 

(31) and high aspect ratios. The correlation coefficients are 
presented in table I. The three correlation curves for ratio of 
axial chord to mean diameter are plotted in figure 3. Using 
one-third of axial chord length for spacing between blade rows, 

(32) the slope for a linear meanline between stage 1 exit and stage 
n exit is 




me, ex 


+ vl 


Pex = Pin 



\ 


(33) tan 0 = Dm, In) _ 3 (Dm, ex Dmjn) ^ 

2(4/3 )E CxJ 16 Ec,.,- 

2 2 

(34) For two-segment meanlines, the lower or upper limit for the 
summation is changed as appropriate for the linear section, 
and the slope for the constant section is set equal to zero. 
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TABLE L-AXIAL CHORD CORRELATION COEFFICIENTS 
\ C X = k\D m + k 2 Df„ + k 3 + k A D~ l . ] 


Aspect 

ratio 

k t x 10 2 

k 2 X 10 3 

k 3 x 

10 2 

k 4 x 10 3 

m -1 

ft " 1 

m 

ft 

m 2 

ft 2 

High 

0.5716 

7.8984 

2.4074 

1.4911 

4.8921 

-0.7498 

-8.0708 

Mid 

3.0740 

9.9323 

3.0274 

1.4714 

4.8275 

-.7411 

-7.9771 

Low 

6.6559 

8.7134 

2.6558 

1.1949 

3.9203 

-.5580 

-6.0063 


Exit annulus area, radius ratio, and hub and tip diameters 
are obtained as follows: 


S 1 

^ i 

II 

5 

Q. 

(41) 

RT ex 


— ]/ cos 8 

,ex v me, ex u n 

(42) 

W 


^ an, ex ” T . 

(43) 


Pex^x.ex 



1 - 

7tZ)2 

^ m, ex; 

1 + 

A 

an, ex 


itDi 


(44) 


&t,ex 



(45) 



(46) 



0 .5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 

Mean diameter, ft 


Absolute and relative flow angles, which are the same for 
each stage, are computed from the last-stage velocities: 


a 2 — tan 


\ K.2./1 


(47) 


Figure 3.— Axial chord correlation. 

Absolute and relative Mach numbers are computed for the last 
stage, where velocities are highest and temperatures are lowest, 
thus making the Mach numbers most severe: 


Wh,L/i — Kt.Ur 


U„ 


(48) 


"2 ,n 


= & 


u,2,n 


+ vi 


(52) 


K.2,„ = K.2 ,i - U n 

(49) 

ft = tan -1 Wu - Un 
1/ 

(50) 

* me, n 


ft = tan" 1 WuXn 
v 

r mp.n 

(51) 


me,n 


w 2 = V< 2 ,, + fL.,< 

(53) 

Wij, = v<,,„ + VL.„ 

(54) 

V,,, = VKL,„ + Vle.n 

(55) 
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iA — wi 

rpf 2,n 2 ,n 

*2 ,n~ r - 1 ex ' T 

2gJc p 

(56) 

V o 

y ^ “> 2 a V 

* u,2,a ayt ^”u,a 

af„ 

(69) 

T T - Kn 

l \,n - L \,n - , 

2 gJc p 

(57) 

Vme,2,a = ~~ &K,a COt a, 

af„ 

(70) 

M V 

abs - 2 s/jgRT ex 

(58) 

V 2 , a = 'Jv 2 uX a+Vle, 2, a 

(71) 

M 

M tt\,2 — / 

^ygRT ex 

(59) 

v\ a 

rji rpif 

I 2,a ~ 1 2,a _ . 

2gJc p 

(72) 

W Xn 

M rel * = . ===== 
' Vt gRTx, n 

(60) 

/ 7* \yKy~~ l ) 
, / 1 2,a\ 

P2,a P2,a 1 rp, 1 
\ l 2,a/ 

(73) 

** ^1 ,n 

M abs,l ““ i 

vygRT\, n 

(61) 

Kx, 2 ,a ~~ Vme,2,a ^05 

(74) 


Efficiency and continuity computations are based solely on 
the meanline flow parameters. Although radial variations of 
velocity in a turbine need not be of the free-vortex type, the 
free-vortex hub and tip velocity-diagram parameters can give 
a qualitative indication of the possible severity of end wall flow 
conditions. With blade speed directly proportional to diameter 

UocD ( 62 ) 

and applying the free-vortex conditions for radial variations 
of the velocity components 


P2,a “ 


P2,a 

RT 2 ,a 


“ an, 2, a 


W 

P2 r a^x,2,a 



1 - 


Aqn,2,a 


1 + 


Aan,2,a 

irDm.in 


(75) 

(76) 


(77) 


and 


V„cc- 

D 


(63) 



(78) 


V me ~ Constant (64) 

the hub and tip velocity diagrams are determined from 
equations (47) to (61). 

In order to establish the flow annulus geometry near the 
turbine inlet, a flow analysis is made at the first-stage exit as 
follows: 


Dhxa = D, Xa (79) 

When a constant annulus is assumed for the first stage, the 
first-stage exit dimensions become the inlet dimensions for that 
turbine. 


II 

¥ 

(65) 

ryt _ ryt 

l 2 ,a — * in 

(66) 

l 

1 

II 

e 

(67) 

A K.a = m 

(68) 


Alternative Calculation Procedure 

The basic calculation procedure described in the previous 
section requires as inputs the inlet and exit mean diameters, 
stator exit angle, and shaft power. Alternatively, the hub or 
tip diameters could be specified as input, and the mean 
diameters computed; the exit radius ratio could be specified 
as input, and the stator exit angle computed; and/or the turbine 
pressure ratio could be specified as input, and the shaft power 
computed. These alternative input options require iterative 
calculation procedures such as described in this section. 
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With hub or tip diameters rather than mean diameters 
specified at the inlet and exit, it is necessary to assume initial 
values for the inlet and exit radius ratios. Initial values for 
inlet and exit mean diameters are then obtained as 



(80) 


when hub diameter is input and 



(81) 


when tip diameter is input. The computation then proceeds 
from equation (1) to equation (44), and the computed exit 
radius ratio is compared with the assumed value. If they are 
not the same (within a given tolerance), then the computed 
value of exit radius ratio is used to calculate a new value for 
exit mean diameter (from eq. (80) or (81)), and the 
computation procedure is repeated until convergence is 
obtained. Then, computation proceeds through equation (77), 
and the computed inlet radius ratio is compared with the 
assumed value. If they are not within tolerance, the computed 
value of inlet radius ratio is used to calculate a new value for 
inlet mean diameter, and the computation procedure is repeated 
from equation (1). This entire procedure is repeated until both 
inlet and exit radius ratios converge in the same calculation 
pass. When hub or tip diameter rather than mean diameter 
is used as input, a two-segment meanline cannot be specified. 

With exit radius ratio rather than stator exit angle specified, 
a value of stator exit angle is assumed for the evaluation of 
equation (17). The computation proceeds through equation 
(41). Equation (44) is then used to compute the exit annulus 
area from the input value of radius ratio, and the exit axial 
velocity is then obtained from equation (43). The density used 
in equation (43), however,' is not consistent with the exit area, 
and equations (43), (42), (33), (36), and (41) must be iterated 
until convergence is obtained. Then, the stator exit angle is 
computed as 




— tan 



(82) 


and compared with the assumed value. If they are not within 
tolerance, the computed value of stator exit angle is used for 


the evaluation of equation (17), and the computation procedure 
just given is repeated until stator exit angle converges. The 
remainder of the computation is then completed. 

With the turbine inlet-total- to exit-static-pressure ratio rather 
than shaft power specified as input, an initial value of turbine 
static efficiency is assumed. Turbine work is then computed 
from 



(83) 


instead of from equation (4). The computation then proceeds 
through equation (38) to compute a static efficiency. If the 
computed and the assumed values are not within tolerance, 
a new value of static efficiency is assumed, and the computation 
is repeated until convergence. The remainder of the com- 
putation is then completed. 


Loss-Coefficient Evaluation 

Ten turbines that were designed for airbreathing engine 
applications and performance tested in component facilities 
were selected to serve as the basis for the loss coefficient 
evaluation. These turbines represent the output of several 
different design systems and cover a wide range of design 
characteristics. Six of the turbines were designed and tested 
by NASA, two were designed and tested by GE Aircraft 
Engines Co. , and two were designed by Pratt & Whitney, one 
of which was tested by Pratt & Whitney and the other by 
NASA. The characteristics of these turbines are presented in 
table II. The number of stages varies from 1 to 5, stage work 
factors (^) from 1.2 to 4.7, stator angles from 51° to 80°, 
exit mean diameters, except for one, from 47 to 69 cm (18 
to 27 in.), and radius ratios from 0.50 to 0.86. All the listed 
total efficiencies except for two were as measured. In those 
two cases, corrections were made to account for the effects 
of cooling air (ref. 4) and a large discrepancy in flow rate from 
design intent (ref. 5). 

These turbines were modeled as closely as possible for the 
subject computer program, and design performance was 
computed over a range of loss coefficients. The variation of 
efficiency with loss coefficient was a direct proportion between 
the quantity (1 — rj) and loss coefficient, K. Calculated 
efficiencies are plotted against measured efficiencies in figure 4 
for a constant loss coefficient of 0.3. As seen, for all cases 
except one, the calculated efficiency was within one point of 
the measured efficiency. In view of the simplistic modeling 
used for this computer program, this agreement is admittedly 
fortuitous. It does illustrate, however, that this program is 
capable of predicting reasonable efficiencies over a wide range 
of airbreathing engine turbine conditions. A loss coefficient 
value of 0.3 is recommended for use in the absence of 
additional information. 
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TABLE II. —TURBINES USED FOR EVALUATION 


Number 

of 

stages 

Average 

stage 

work 

factor, 

* 

Average 

stator 

angle, 

deg 

Exit 

mean 

diameter 

Exit 

radius 

ratio 

Measured 

efficiency 

Ref- 

erence 

cm 

in. 

1 

1.7 

67.0 

66.0 

26.0 

0.73 

0.923 

6 

2 

1.7 

64.5 

66.0 

26.0 

.63 

.932 

7 

1 

1.2 

63.4 

20.3 

8.0 

.68 

.917 

8 

1 

1.9 

73.1 

47.0 

18.5 

.85 

.886 

9 

3.5 

4.0 

51.4 

49.8 

19.6 

.50 

.880 

5 

4.5 

4.7 

66.6 

49.8 

19.6 

.57 

.855 

10 

i 

1.7 

75.0 

46.7 

18.4 

.86 

.895 

11 

1 

1.6 

79.6 

56.4 

22.2 

.85 

.901 

12 

2 

1.3 

71.6 

69.3 

27.3 

.82 

.898 

4 

5 

2.6 

61.0 

63.5 

25.0 

.62 

.920 ! 
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Measured efficiency 


Figure 4.— Comparison of measured and calculated total efficiencies. Loss 
coefficient, 0.3. 


Description of Input and Output 

This section presents a detailed description of the program 
input, normal output, and error messages. Included in the input 
and output sections are several example cases illustrating the 
use of the program and the various options. 

Input 

The program input, a sample of which is presented in table HI, 
consists of a title record and the required physical data and 
option indicators in namelist form. The title, which is printed 
as a heading on the output listing, can contain up to 77 
characters located anywhere in columns 2 to 78 on the title 
record. A title, even if it is left blank, must be the first record 
of the input data. Additional titles can be used to identify 
different cases being run in the same data file. This is done 
by placing a title in front of the data for the particular case 
and using the option indicator itit as subsequently described . 

The physical data and option indicators are input in data 
records having the namelist name input. The variables and 
indicators that compose input and the proper units are as 
follows. These must be input for all cases except where 
otherwise indicated. Either the SI units or the U.S. customary 
units shown below may be used. 

PTiN inlet total pressure, N/cm 2 ; lb/in. 2 
TTIN inlet total temperature, K; °R 
mu gas viscosity, N* sec/m 2 ; lb/sec -ft 
r gas constant, J/kgvK; ft • Ibf/lbm* °R 

gam specific heat ratio 

din inlet diameter-hub or mean or tip value as 
specified by the indicator idiam, cm; in. 

dex exit diameter— hub or mean or tip value as 

specified by the indicator idiam, cm; in. 

rrex exit radius ratio; rrex may be omitted in the 

case where both idiam — 2 and ialph = 0; rrex 
is used as first trial value when ialph = 0 and 
idiam = 1 or 3 

rpm rotative speed, rad/sec; rpm 

pow shaft power— omit when ipr = 1, kW; hp 

w mass flow rate, kg/sec; lb/sec 


TABLE III. -SAMPLE INPUT 

THESE ARE SAMPLE CASES FOR THIS PROGRAM. U.S. CUSTOMARY UNIT ARE USED , 
&INPUT PTIN=113.1,TTIN=266Q. , MU= . 376E- 04 ,R=53 . 37 , GAM^l .302,DIN=22. ,DEX=24. , 
RPM=11400. ,P0W=12900. ,W=53 . 5 , ALPHA=65 . ,KL0SS= . 35 , NMIN=2, NMAX=2 , E=1 .2, 
IALPH-O , IDIAM=2 , I VD= I , IEV=0 , IPR=0,IU=2 SEND 
& INPUT IDIAM=1,RREX=.8,DIN=20. ,DEX=21. , I'T'IT-1 SEND 

THIS IS AN ADDITIONAL TITLE CALLED FOR BY ITIT IN PREVIOUS CASE DATA 
&INPUT IALPH” I , IEV=1 , SEND 
SINPUT IALPH=0 / IEV=0 , IVD=2 SEND 
&INPUT IVD=3 SEND 
SINPUT IPR=1 , IVD-l , PRTS-3 . 66 SEND 
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alpha stator exit angle from axial direction; alpha is 
used as first trial value when ialph = 1 , deg 

alphao turbine inlet flow angle; required only when 
kalpho ■= 2, deg 

vuidvu ratio of rotor inlet swirl to total change in 
swirl; input only when ivd = 5 

KLOSS turbine loss coefficient; a value of 0.3 is 
recommended in the absence of additional 
information 

nmin minimum number of stages for which the 
calculations are performed 

nmax maximum number of stages for which the 

calculations are performed; results are obtained 
for all stage numbers between nmin and nmax 

nmid stage number at which meanline changes slope; 
may be omitted when imid — 0 

e squared ratio of stage-exit to stage-average 

meridional velocities 

prts turbine inlet-total to exit-static-pressure ratio; 
omit when ipr — 0 

ialph indicates whether stator exit angle or turbine 
exit radius ratio is specified: 
ialph = 0— turbine is designed for specified 
alpha 

ialph = 1— turbine is designed for specified 

RREX 

idiam indicates whether input diameters are hub, mean, 
or tip values: 

IDIAM = 1— input diameters are hub values 
idiam == 2— input diameters are mean values 
IDIAM — 3— input diameters are tip values 

ivd indicates type of velocity diagram used: 
ivd = 1 — symmetrical diagrams 
ivd = 2— zero exit swirl diagrams 
ivd = 3— impulse diagrams 

ivd = 4— zero exit swirl diagrams if < 2.0 
and impulse diagrams if $ > 2.0 
ivd = 5— ratio of rotor inlet swirl to total 
change in swirl is input as vuidvu 

itit indicates use of titles in addition to that required 
as first line of data package: 

itit = 1— title line precedes next data set; 
must be input for each additional title because 
itit is automatically restored to zero after 
each title is read 

iev indicates use of exit vanes: 
iev = 0— no exit vanes 

iev = 1— exit vanes are used to turn turbine 
exit flow to axial direction 


ipr indicates whether shaft power or pressure ratio 
is specified: 

ipr = 0— shaft power is input 

ipr — 1 —turbine inlet-total- to exit-static- 
pressure ratio is input 

iu indicates type of units used for input and output: 

iu = 1— SI units 
iu — 2— U.S. customary units 
kalpho indicates turbine-inlet flow angle option: 

kalpho — 0— turbine-inlet flow is axial (default) 

kalpho = 1— turbine-inlet flow angle equals 
stage-exit flow angle 

kalpho = 2— turbine-inlet flow angle is input 
as alphao 

iar indicates blading aspect ratio: 

iar ■== 1— high aspect-ratio blading 
iar = 2— mid aspect-ratio blading (default) 
iar = 3— low aspect-ratio blading 
imid indicates meanline shape: 

imid = 0— meanline linear from stage 1 to 
stage N (default) 

imid — 1— meanline constant from stage 1 to 
stage nmid and then linear to stage N 

imid = 2— meanline linear from stage 1 to 
stage nmid and then constant to stage N 

The first line of the input file shown in table III is the 
mandatory title card, which can contain any desired message. 
The next three lines are the first data set, which contains all 
required inputs. This first case represents computation in 
accordance with the basic calculation procedure described 
previously. Data input for subsequent cases need only include 
those values that differ from previous ease data. The fifth line 
is the second data set and represents the option where hub 
diameter is input. Also, the second case data specify that a 
title, which is the sixth data line, precede the third case data. 
Cards 7 to 10 represent four additional cases illustrating the 
use of different input options. The output corresponding to 
this sample input is described in the following section. 

Output 

The program output consists of title headings, the input 
variables, and computed results. This section presents normal 
output. Error message output is described in the next section. 

Table IV presents the output that corresponds to the sample 
input shown in table III. The top line of output is a program 
identification title that is automatically printed. The second 
line is the title card message. The next four lines are the input 
variables and their associated values for the first data set. The 
input variable names are spelled out. The units for the input 
variable values are as described in the “Input” section. The zeros 
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printed under EXIT RADIUS RATIO and T-S PRESSURE 
RATIO indicate that these are computed for this case and not 
specified by the input. The fact that the input diameters are 
mean diameters is indicated by the MN in the variable name. 
Hub and tip diameters would be indicated by HB and TP, 
respectively. These four lines of output are printed for each 
new data set. 

The next group of eight lines is the computation results 
satisfying the input requirements. The output parameters are 
spelled out and are self-explanatory. These temperatures, 
pressures, velocities, and angles are meanline values. On the 
first line are the number of stages, the stage work factor, the 
Reynolds number, and the diagram type, which is symmetrical 
for this first case. The remainder of the output includes exit 
and inlet tip and hub diameters in the first column, exit total 
and static temperatures and pressures in the second column, 
total and static efficiencies and velocity diagram angles in the 
third column, and first- and last-stage blade speeds and last- 
stage absolute velocity components in the last column. The 
last line of this output group presents the last-stage absolute 
and relative Mach numbers. 

The next group of four lines is the hub and tip free-vortex 
values of Mach numbers and angles. The four columns present 
the stator-exit absolute, rotor-inlet relative, rotor-exit relative, 
and stage-exit absolute values. The last line of output for each 
case is the meanline slope based on the specified (mid in this 
case) aspect-ratio blading. 


After the computations for each input case are completed, 
the input data for the next case are printed. The second case 
presented here is that where the specified diameters are hub 
values. The third input case in table IV is preceded by an 
additional title line. This causes the next output to begin at 
the top of a new page with the program identification and the 
title message. 

The third case is for a specified exit radius ratio rather than 
for a specified stator exit angle. This is indicated by a zero 
appearing under stator ex ang in the row of input variables. 
Exit vanes are included in this case, as indicated by with exit 
vanes printed after number of stages on the first line of result 
output. The fourth and fifth cases, as indicated by the top line 
of the fourth column of result. output, are for zero exit swirl 
and impulse diagrams, respectively. Turbine pressure ratio 
rather than shaft power is specified in the sixth case. For this 
case, the input pressure ratio was chosen to be the same as 
the computed value obtained for the second case. As should 
be expected, the program converges to identical solutions for 
both cases. An extra line of output consisting of shaft power 
and specific work is printed for the case where pressure ratio 
is specified as input. 

Error Messages 

The program contains seven output messages. Five of these 
indicate the nonexistence of a meanline solution satisfying 


table iv.-sample output 


TURBINE VELOCITY DIAGRAM ANALYSIS 

THESE ARE SAMPLE CASES FOR THIS PROGRAM. U.S. CUSTOMARY UNIT ARE USED. 


SHAFT 

POWER 

12900.0 


MASS 

INLET 

INLET 

ROTATIVE 

INLET 

EXIT 

EXIT 

STATOR 

GAS 

FLOW 

TEMP 

PRESS 

SPEED 

MN DIA 

MN DIA 

RADIUS 

RATIO 

EX ANG 

CONST 

53.50 

2660.00 

113.10 

11400.00 

22.00 

24.00 

0.0000 

65.00 

53.37 


STAGES* 2 

EXIT TIP DIAMETER = 27.18 
EXIT HUB DIAMETER = 20.82 
EXIT RADIUS RATIO =0.7658 
INLET TIP DIAMETER* 24.16 
INLET HUB DIAMETER* 19.84 
INLET RADIUS RATI0=0.8213 
LAST STG Ml ABS =0.7810 


STAGE WORK FACTOR* 1.63 
EXIT TOTAL TEMP =2083.64 
EXIT STATIC TEMP =2030.85 
EXIT TOTAL PRESS = 34.85 

EXIT STATIC PRESS= 31.20 
T-T PRESS RATIO = 3.245 

T-S PRESS RATIO = 3.625 

LAST STG Ml REL = 0.3708 


REYNOLDS NO . =0 . 1552E+07 
STATOR EXIT ANGLE= 65.00 
STAGE EXIT ANGLE =-27.10 
ROTOR INLET ANGLE* 27.10 
ROTOR EXIT ANGLE =-65.00 
TOTAL EFFICIENCY =0.907 
STATIC EFFICIENCY* 0 .839 
LAST STG M2 REL *0.8263 


HUB: 

TIP: 


LAST STG Ml ABS =0.8915 
STATOR EXIT ANGLE = 67.98 
LAST STG Ml ABS =0.7010 
STATOR EXIT ANGLE = 62.16 


LAST STG Ml REL = 0.4862 
ROTOR INLET ANGLE * 46.57 
LAST STG Ml REL = 0.3277 
ROTOR INLET ANGLE = 2.51 


LAST STG M2 REL =0.7684 
ROTOR EXIT ANGLE =-63.51 
LAST * STG M2 REL =0.8590 
ROTOR EXIT ANGLE =-66.51 


STAGE 1- 2 MEANLINE SLOPE = 14.85 DEG BASED ON MID ASPECT-RATIO BLADING 


HEAT 

GAS 

TURBINE 

AXIAL 

T-S 

CAPAC 

VISCOSITY 

LOSS 

VEL SQ 

PRES 

RATIO 


C0EF 

RATIO 

RATI 

1.302 

0 . 376E-04 

0.350 

1.200 

0 . 


DIAGRAMS ARE SYMMETRICAL 
FIRST STAGE MEAN SPEED* 1094.32 
LAST STAGE MEAN SPEED = 1193.81 
LAST STAGE INLET SWIRL* 1567.87 
LAST STAGE EXIT SWIRL = -374.06 
LAST STAGE MERID VEL0C= 800.93 
EXIT MERID MACH NUMBER* 0.3759 
LAST STG M2 ABS * 0.4149 

LAST STG M2 ABS = 0.3980 

STAGE EXIT ANGLE =-30.53 

LAST STG M2 ABS = 0.3757 

STAGE EXIT ANGLE *-24.31 


SHAFT 

POWER 

12900.0 


MASS 

INLET 

INLET 

ROTATIVE 

INLET 

EXIT 

FLOW 

TEMP 

PRESS 

SPEED 

HB DIA 

HB DIA 

53, 50 

2660.00 

113.10 

11400.00 

20.00 

21.00 


EXIT STATOR GAS 
RADIUS EX ANG CONST 
RATIO 

0.0000 65.00 53.37 


HEAT 

GAS 

TURBINE 

AXIAL 

T-S 

CAPAC 

VISCOSITY 

LOSS 

VEL sq 

PRES 

RATIO 


COEF 

RATIO 

RATI 

1 .302 

0 . 376E-04 

0.350 

1.200 

0 . 


STAGES* 2 

EXIT TIP DIAMETER = 27.32 
EXIT HUB DIAMETER = 21.00 
EXIT RADIUS RATIO =0.7685 
INLET TIP DIAMETER= 24.29 
INLET HUB DIAMETER* 20.00 
INLET RADIUS RATI0=0.8233 
LAST STG Ml ABS =0.7798 


STAGE WORK FACTOR* 1.61 
EXIT TOTAL TEMP =2083.64 
EXIT STATIC TEMP =2031.50 
EXIT TOTAL PRESS = 34.87 

EXIT STATIC PRESS* 31.27 
T-T PRESS RATIO = 3.243 

T-S PRESS RATIO * 3.617 

LAST STG Ml REL = 0.3682 


REYNOLDS NO . =0 . 1542E+Q7 
STATOR EXIT ANGLE* 65.00 
STAGE EXIT ANGLE =-26.48 
ROTOR INLET ANGLE* 26.48 
ROTOR EXIT ANGLE *-65.00 
TOTAL EFFICIENCY =0.907 
STATIC EFFICIENCY*0 .840 
LAST STG M2 REL =0.8249 


DIAGRAMS ARE SYMMETRICAL 
FIRST STAGE MEAN SPEED* 1101.56 
LAST STAGE MEAN SPEED = 1201.87 
LAST STAGE INLET SWIRL* 1565.51 
LAST STAGE EXIT SWIRL = -363.64 
LAST STAGE MERID VEL0C* 799.73 
EXIT MERID MACH NUMBER* 0.3753 
LAST STG M2 ABS = 0.4123 


HUB: LAST STG Ml ABS =0.8883 
STATOR EXIT ANGLE = 67.94 
TIP: LAST STG Ml ABS =0.7007 
STATOR EXIT ANGLE = 62.20 


LAST STG Ml REL = 
ROTOR INLET ANGLE = 
LAST STG Ml REL * 
ROTOR INLET ANGLE = 


0.4805 

46.03 

0.3271 

1.97 


LAST STG M2 REL =0.7664 
ROTOR EXIT ANGLE =-63.48 
LAST STG M2 REL =0.8580 
ROTOR EXIT ANGLE =-66.52 


STAGE 1- 2 MEANLINE SLOPE = 14.89 DEG BASED ON MID ASPECT-RATIO BLADING 


LAST STG M2 ABS 
STAGE EXIT ANGLE 
LAST STG M2 ABS 
STAGE EXIT ANGLE 


= 0.3944 
=-29.82 
= 0.3735 
=-23.77 
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TABLE IV.— Concluded. 


TURBINE VELOCITY DIAGRAM ANALYSIS 


THIS IS AN ADDITIONAL TITLE CALLED FOR BY ITIT IN PREVIOUS CASE DATA 


SHAFT 

MASS 

INLET 

INLET 

ROTATIVE 

INLET 

EXIT 

EXIT 

STATOR 

GAS 

HEAT 

GAS TURBINE 

AXIAL 

T-S 

POWER 

FLOW 

TEMP 

PRESS 

SPEED 

HB DIA 

HB DIA 

RADIUS 

RATIO 

EX ANG 

CONST 

CAPAC 

RATIO 

VISCOSITY LOSS 

COEF 

VEL SQ 
RATIO 

PRESS 

RATIO 

12900.0 

53.50 

2660.00 

113.10 

11400.00 

20.00 

21.00 

0.8000 

0 . 00 

53.37 

1.302 

0 . 376E-04 0.350 

1.200 

0.000 


STAGES* 2 W/EXT VANES STAGE WORK FACTOR* 1.67 REYNOLDS NO .*0 . 1569E+07 DIAGRAMS ARE SYMMETRICAL 


EXIT TIP DIAMETER = 26.25 EXIT TOTAL TEMP *2083.64 STATOR EXIT ANGLE* 59.54 FIRST STAGE MEAN SPEED* 1082.83 

EXIT HUB DIAMETER * 21.00 EXIT STATIC TEMP =2014.59 STAGE EXIT ANGLE =-23.12 LAST STAGE MEAN SPEED = 1175.15 

EXIT RADIUS RATIO =0.8000 EXIT TOTAL PRESS = 34.80 ROTOR INLET ANGLE* 23.12 LAST STAGE INLET SWIRL* 1569.27 

INLET TIP DIAMETER* 23.54 EXIT STATIC PRESS* 30.10 ROTOR EXIT ANGLE =-59.54 LAST STAGE EXIT SWIRL = -394.12 

INLET HUB DIAMETER* 20.00 T-T PRESS RATIO = 3.250 TOTAL EFFICIENCY =0.906 LAST STAGE MERID VELOC* 1011.01 

INLET RADIUS RATIO=0.8497 T-S PRESS RATIO = 3.758 STATIC EFFICIENCY=0 . 81 9 EXIT MERID MACH NUMBER* 0.4764 

LAST STG Ml ABS =0.8263 LAST STG Ml REL * 0.4555 LAST STG M2 REL =0.8797 LAST STG M2 ABS * 0,5113 

HUB: LAST STG Ml ABS =0.9135 LAST STG Ml REL = 0.5370 LAST STG M2 REL *0.8255 LAST STG M2 ABS * 0,4827 

STATOR EXIT ANGLE = 62.40 ROTOR INLET ANGLE = 37.99 ROTOR EXIT ANGLE =-58.19 STAGE EXIT ANGLE =-25.66 

TIP: LAST STG Ml ABS =0.7603 LAST STG Ml REL * 0.4186 LAST STG M2 REL =0.8945 LAST STG M2 ABS = 0.4656 

STATOR EXIT ANGLE = 56.84 ROTOR INLET ANGLE = 6.59 ROTOR EXIT ANGLE =-60.93 STAGE EXIT ANGLE =-21.02 


STAGE 1- 2 MEANLINE SLOPE = 13.98 DEG BASED ON MID ASPECT-RATIO BLADING 


SHAFT 

MASS 

INLET 

INLET 

ROTATIVE 

INLET 

EXIT 

EXIT 

STATOR 

GAS 

HEAT 

GAS 

TURBINE 

AXIAL 

T-S 

POWER 

FLOW 

TEMP 

PRESS 

SPEED 

HB DIA 

HB DIA 

RADIUS 

RATIO 

EX ANG 

CONST 

CAPAC 

RATIO 

VISCOSITY 

LOSS 

COEF 

VEL SQ 
RATIO 

PRESS 

RATIO 

12900.0 

53.50 

2660.00 

113.10 

11400.00 

20.00 

21.00 

0.0009 

65.00 

53.37 

1.302 

0 . 376E-04 

0.350 

1.200 

0.000 


STAGES* 2 STAGE WORK FACTOR* 1.67 REYNOLDS NO . =0 . 1568E+07 DIAGRAMS ARE ZERO EX SWRL 


EXIT TIP DIAMETER = 26.34 EXIT TOTAL TEMP *2083.64 STATOR EXIT ANGLE* 65.00 FIRST STAGE MEAN SPEED* 1083.42 

EXIT HUB DIAMETER * 21.00 'EXIT STATIC TEMP =2015.77 STAGE EXIT ANGLE * 0.00 LAST STAGE MEAN SPEED = 1177.29 

EXIT RADIUS RATIO =0.7974 EXIT TOTAL PRESS = 34.43 ROTOR INLET ANGLE* 40.62 LAST STAGE INLET SWIRL* 1962.12 

INLET TIP DIAMETER* 23.56 EXIT STATIC PRESS* 29.85 ROTOR EXIT ANGLE =-52.15 LAST STAGE EXIT SWIRL * 0.00 

INLET HUB DIAMETER* 20.00 T-T PRESS RATIO = 3.285 TOTAL EFFICIENCY =0.899 LAST STAGE MERID VELOC* 1002.33 

INLET RADIUS RATI0=0.8488 T-S PRESS RATIO * 3.789 STATIC EFFICIENCY=0 . 815 EXIT MERID MACH NUMBER* 0,4722 

LAST STG Ml ABS *1.0042 LAST STG Ml REL = 0.5592 LAST STG M2 REL *0.7284 LAST STG M2 ABS * 0.4722 

HUB: LAST STG Ml ABS =1.1294 LAST STG Ml REL = 0,6998 LAST STG M2 REL =0.6523 LAST STG M2 ABS = 0.4298 

STATOR EXIT ANGLE = 67.52 ROTOR INLET ANGLE * 51.90 ROTOR EXIT ANGLE =-48.78 STAGE EXIT ANGLE = 0.00 

TIP: LAST STG Ml ABS =0,9106 LAST STG Ml REL * 0.4681 LAST STG M2 REL =0.7507 LAST STG M2 ABS = 0.4298 

STATOR EXIT ANGLE * 62.58 ROTOR INLET ANGLE = 26.36 ROTOR EXIT ANGLE =-55.07 STAGE EXIT ANGLE = 0.00 


STAGE 1- 2 MEANLINE SLOPE = 14.19 DEG BASED ON MID ASPECT-RATIO BLADING 


SHAFT 

MASS 

INLET 

INLET 

ROTATIVE 

INLET 

EXIT 

EXIT 

STATOR 

GAS 

HEAT 

GAS 

TURBINE 

AXIAL 

T-S 

POWER 

FLOW 

TEMP 

PRESS 

SPEED 

HB DIA 

HB DIA 

RADIUS 

RATIO 

EX ANG 

CONST 

CAPAC 

RATIO 

VISCOSITY 

LOSS 

COEF 

VEL SQ 
RATIO 

PRESS 

RATIO 

12900,0 

53.50 

2660.00 

113.10 

11400.00 

20.00 

21.00 

0 . 0000 

65.00 

53.37 

1.302 

0 . 376E-04 

0.350 

1.200 

0.000 


STAGES* 2 STAGE WORK FACTOR* 1.68 REYNOLDS NO. =0.I576E+07 DIAGRAMS ARE IMPULSE 

EXIT TIP DIAMETER = 26.11 EXIT TOTAL TEMP =2083.64 STATOR EXIT ANGLE* 65.00 FIRST STAGE MEAN SPEED* 1078.13 

EXIT HUB DIAMETER * 21.00 EXIT STATIC TEMP =1999.26 STAGE EXIT ANGLE * 10,46 LAST STAGE MEAN SPEED = 1171.57 

EXIT RADIUS RATIO =0.8044 EXIT TOTAL PRESS * 33.80 ROTOR INLET ANGLE* 44.42 LAST STAGE INLET SWIRL* 2157,44 

INLET TIP DIAMETER* 23.35 EXIT STATIC PRESS* 28.28 ROTOR FXTT ANGLE =-44 . 42 LAST STAGE EXIT SWIRL = 185.70 

INLET HUB DIAMETER* 20.00 T-T PRESS RATIO * 3.346 TOTAL EFFICIENCY ■0.887 LAST STAGE MERID VELOC* 1102.11 

INLET RADIUS RATIO=0.8566 T-S PRESS RATIO * 3.999 STATIC EFFICIENCY^, 788 EXIT MERID MACH NUMBER* 0.5213 

LAST STG Ml ABS =1.1222 LAST STG Ml REL = 0,6640 LAST STG M2 REL =0.6995 LAST STG M2 ABS = 0.5287 

HUB: LAST STG Ml ABS =1.2610 LAST STG Ml REL * 0.8199 LAST STG M2 REL *0.6168 LAST STG M2 ABS = 0.4844 

STATOR EXIT ANGLE * 67,42 ROTOR INLET ANGLE = 53.81 ROTOR EXIT ANGLE =-39.74 STAGE EXIT ANGLE = 11.70 

TIP: LAST STG Ml ABS =1.0182 LAST STG Ml REL = 0.5561 LAST STG M2 REL =0.7135 LAST STG M2 ABS * 0.4808 

STATOR EXIT ANGLE = 62.67 ROTOR INLET ANGLE * 32.78 ROTOR EXIT ANGLE =-48.35 STAGE EXIT ANGLE = 9.45 

STAGE 1- 2 MEANLINE SLOPE = 14.18 DEG BASED ON MID ASPECT-RATIO BLADING 


SHAFT 

MASS 

INLET 

INLET 

ROTATIVE 

INLET 

EXIT 

EXIT 

STATOR 

GAS 

HEAT 

GAS 

TURBINE 

AXIAL 

T-S 

POWER 

FLOW 

TEMP 

PRESS 

SPEED 

HB DIA 

HB DIA 

RADIUS 

RATIO 

EX ANG 

CONST 

CAPAC 

RATIO 

VISCOSITY 

LOSS 

COEF 

VEL SQ 
RATIO 

PRESS 

RATIO 

0.0 

53,50 

2660.09 

113.10 

11400.00 

20.00 

21 . 00 

0,0000 

65,00 

53.37 

1.302 

0 . 376E-04 

0.350 

1.200 

3.660 


STAGES* 2 STAGE WORK FACTOR* 1.62 REYNOLDS NO . =0 . 1542E+07 DIAGRAMS ARE SYMMETRICAL 


EXIT TIP DIAMETER = 27.35 EXIT TOTAL TEMP =2079.52 STATOR EXIT ANGLE* 65.00 FIRST STAGE MEAN SPEED* 1101.57 

EXIT HUB DIAMETER = 21.00 EXIT STATIC TEMP =2026.66 STAGE EXIT ANGLE =-26.78 LAST STAGE MEAN SPEED = 1202.59 

EXIT RADIUS RATIO =0.7677 EXIT TOTAL PRESS = 34.53 ROTOR INLET ANGLE = 26.78 LAST STAGE INLET SWIRL* 1572.70 

INLET TIP DIAMETER* 24.29 EXIT STATIC PRESS* 30.90 ROTOR EXIT ANGLE =-65.00 LAST STAGE EXIT SWIRL = -370.11 

INLET HUB DIAMETER* 20.00 T-T PRESS RATIO = 3.276 TOTAL EFFICIENCY =0.907 LAST STAGE MERID VELOC* 803.40 

INLET RADIUS RATIO=0.8233 T-S PRESS RATIO = 3.660 STATIC EFFICIENCY*0 .840 EXIT MERID MACH NUMBER* 0.3775 

SHAFT POWER = 12992.2 SPECIFIC WORK * 171.68 

LAST STG Ml ABS =0.7840 LAST STG Ml REL = 0.3711 LAST STG M2 REL =0.8297 LAST STG M2 ABS * 0.4156 

HUB: LAST STG Ml ABS *0.8937 LAST STG Ml REL = 0.4852 LAST STG M2 REL =0.7713 LAST STG M2 ABS * 0.3981 

STATOR EXIT ANGLE = 67.95 ROTOR INLET ANGLE = 46.25 ROTOR EXIT ANGLE =-63.50 STAGE EXIT ANGLE =-30.16 

TIP: LAST STG Ml ABS =0.7042 LAST STG Ml REL = 0.3289 LAST STG M2 REL =0.8626 LAST STG M2 ABS = 0.3764 

STATOR EXIT ANGLE = 62.18 ROTOR INLET ANGLE = 2.30 ROTOR EXIT ANGLE *-66.51 STAGE EXIT ANGLE *-24.04 


STAGE 1- 2 MEANLINE SLOPE = 14.99 DEG BASED ON MID ASPECT-RATIO BLADING 

the specified input requirements, while the other two indicate It indicates that the turbine specific work requirement is greater 

the nonexistence of a hub or tip solution satisfying free- vortex than the energy available in the gas. Therefore, either the 

conditions. These messages are presented in this section, and specified shaft power must be decreased or the specified flow 
their causes are discussed. must be increased. 

(1) insufficient energy— This message is caused by the (2) insufficient ideal energy— This message is caused 

computed turbine exit total temperature being less than zero. by the computed ideal energy being more than that available 
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from an infinite expansion of the gas. It indicates that the 
computed efficiency is too low to yield a valid solution. 
Corrective action includes decreasing power, increasing flow, 
using more stages, or perhaps using a different velocity 
diagram. 

(3) negative tex— This message is caused by the computed 
turbine exit static temperature being less than zero. It indicates 
a low value of turbine exit total temperature and/or a high value 
of turbine exit velocity. Corrective action could be the same 
as for message (2) or decreasing a too high value of stage exit 
to average meridional kinetic energy ratio (E). 

(4) insufficient exit area— This message is caused by 
the computed exit area being larger than that available in the 
turbine. Such a situation can be remedied in many ways, 
including increasing exit diameter, decreasing stator exit angle, 
increasing inlet pressure, decreasing mass flow, and increasing 
stage exit to average meridional kinetic energy ratio. 

(5) insufficient inlet area— This message is caused by 
the computed inlet area being larger than that available in the 
turbine. Corrective measures are the same as for message (4). 

(6) negative t at stator hub exit— This message is 
caused by the stator-exit hub static temperature being less than 
zero as a result of a very high absolute velocity. Although a 
non-free-vortex design (not available from this program) might 
not result in a similar condition, it is extreme enough to indicate 
that the mean-section design may require change such as a 
larger diameter, smaller stator angle, or more stages. This 
condition does not interrupt the mean section computations. 

(7) negative t at rotor tip exit— This message is caused 
by the rotor-exit tip static temperature being less than zero 
as a result of a very high absolute velocity. Corrective 
measures are the same as for message (6) above. 

National Aeronautics and Space Administration 
Lewis Research Center, Cleveland, Ohio 44135 
August 12, 1991 
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